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Abstract 



Process of preheating for supersymmetric hybrid model of inflation is generally analyzed 

in two different ways known as parametric resonance and tachyonic preheating. In a common 

frame-work, we study both the processes from the decay or growth of homogeneous inflaton 

field and inhomogeneous waterfall field. We find that these two processes in SUSY hybrid 

F-term inflation are not parameter independent; rather one process will be more preferable 

than the other depending on the values of parameters and Fourier mode of the waterfall field. 

Parameters of the inflationary potential are constrained from the observed CMB data and 

! these constrains help us to identify the process of preheating responsible for this model. 

Ph. 



1 Introduction 



Theory of inflation has been successful so far in describing the origin of observed cosmic microwave 
background (CMB) anisotropics. Since last decade COBE and WMAP has reached a certain level 
^ \ in measuring the values of the cosmological parameters. Planck experiment is aimed at increasing 
the accuracy of those data in recent future. As the accuracy of measurement increased some models 
of inflation have been ruled out, and some models became more successful. 

Standard hybrid irmation[I] has been ruled out just after the release of WMAP first year data[2]. 
q ! But supersymmetric version of hybrid inflation [5] survived because it gives a red spectrum as 
<N observed by WMAP. At the same time this model suggests a negligible amplitude of gravity wave 
• • ' which is to be verified in upcoming Planck experiment jl]. 

.£h ' Preheating is a process of rapid decay of inflaton's energy to other forms just after end of 
^ . inflation. There are different mechanisms of preheating. Firstly there is the old and standard 
c5 '■ concept of preheating by broad parametric resonance [5]. Possibility of this mechanism to be 
efficient has been studied in the context of susy- hybrid model [6] . Another process of preheating is 
called tachyonic preheating, which became more important than this standard one in susy-hybrid 
models because of the negative curvature of the potential after the end of inflation. 

Analysis of preheating in susy hybrid model has been an important field of study since last one 
decade. For preheating process in hybrid inflation [5], the trajectory of the fields were analyzed 
taking the inflaton field and the waterfall field as homogeneous. Only small inhomogeneous pertur- 
bation was considered in the background of a large homogeneous component of the waterfall field. 
Parametric resonance was analyzed, taking fully inhomogeneous waterfall field in ref.[7]. When 
tachyonic preheating is analyzed, waterfall field is generally taken to be completely inhomogeneous. 
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In ref.[6j for SUSY hybrid inflation an effectively one field potential was derived by assuming a 
trajectory of the homogeneous fields and neglecting the oscillation of inflaton. That form of po- 
tential was used in ref.[8] to analyze the tachyonic preheating process in susy-hybrid inflation case. 
Important numerical simulations have been done with LATTICEEASYJ5] and DEFROST [TO] code. 
These codes helped greatly to study the growth of the perturbations with time. But often hard to 
identify the process responsible for this growth. 

In this work we did not assume any effective form of the susy-hybrid potential. Rather we have 
taken the full two-field potential. Inflaton is considered completely homogeneous, but waterfall field 
is considered fully inhomogeneous, i.e. no background homogeneous component is assumed. So we 
do observe an effect of oscillating potential coming from the oscillation of inflaton field. Also, the 
negative curvature part in the potential boosts tachyonic growth of the modes of waterfall field. We 
will show here that in certain range of parameter space this oscillating potential will lead to efficient 
parametric resonance for some Fourier modes of waterfall field, while others modes may experience 
tachyonic growth. So, tachyonic preheating and parametric resonance are not two different processes 
to happen at two different times; rather at the same time different Fourier modes of waterfall field 
may encounter different processes. But for other ranges of parameters where amplitude of inflaton's 
oscillation decays down very quickly, parametric resonance is not so efficient. Therefore tachyonic 
preheating will be the main dominating process for most of the modes in that parameter range. 

So, firstly from the results of our analysis we can identify two different processes of preheating at 
the same time. Secondly, we do show that preheating process in susy-hybrid inflation is a parameter 
dependent process, not as it was expected in ref. [5J. Depending on the values of parameter the 
process of preheating for different modes of the waterfall field gets changed. 

We will be focusing on the F-term susy-hybrid inflation in our work. Inflation occurs due to 
the 1-loop correction in the potential along certain flat direction. In every model of inflation some 
amount of parameter tuning is necessary to fit the observed data. So here also certain values of 
the parameters are favored for latest CMB observation. Can those tuned parameter say anything 
about the mechanism of preheating in susy-hybrid model? In this report we will try to address this 
question. 

In first section we will first describe the SUSY hybrid inflation formalism and show how to 
calculate the observed CMB quantities out of it. After that different mechanisms of preheating will 
be discussed for this model. First we will describe the process of parametric resonance with fully 
inhomogeneous waterfall field. We will show that how different modes of the waterfall filed behave 
differently under parametric resonance. Condition for efficient parametric resonance will also be 
discussed in this section. Next we will briefly outline the tachyonic preheating process. It is mostly 
impossible to workout the fields' amplitude and occupation number analytically. So the analysis 
will be done mostly numerically. 

We present here a common formalism to describe both the processes simultaneously. Growth 
of occupation number of different modes of the waterfall field gives the hints of what process is 
dominating at that situation. There is a technical problem of considering Hubble parameter (H(t)) 
in solving the equation of motion of these fields. Hubble parameter depends on the spatial average of 
the inhomogeneous filed. In mean field approximation equation of motion also involves the spatial 
average of waterfall field. So to solve this exactly one need to solve infinite number of coupled 
equations. This technical problem forces many authors to consider both fields as homogeneous. But 
here to overcome this problem we assume that the spatial average of the square of the inhomogeneous 
waterfall field initially grows exponentially as tachyonic preheating and settles at the minima of the 
potential. Assuming this we see the evolution of Hubble parameter. And then using that we solve 
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the different modes of the waterfall field. After that we will show how different choices of parameters 
of the model can affect the process of preheating. We will also show the mode by mode pattern of 
the growth in occupation number for different parameter choice. We find a clear distinction between 
the patterns of the growth of occupation number for different choice of parameters. 

In the last section we will find out what values of observed CMB quantities can this model 
give to us. To fit with the observed values of these quantities parameters of the model will be 
tuned. Then we show what can those required values of parameters tell us about the mechanism of 
preheating. Choice of values of these parameters depends on the choice of 1-loop correction in the 
effective potential. There are lots of different models in the susy-hybrid class where different forms 
of one loop correction are used. Some of them have observational motivation and some of them 
have theoretical motivation. Some times extra symmetries are also imposed in the superpotential. 
For these different varieties of models tuned parameters will be different. So depending on that 
conclusion may differ from model to model. Here we take the old and standard susy-hybrid 1-loop 
correction just as an example. Although the analysis of preheating process does not depend on the 
1-loop potential. It is because this effective potential is valid only in the time of inflation. After the 
end of inflation a soft susy breaking mass term drives inflaton to the minima, and the value of this 
term is very small compared to the 1-loop correction term. So our analysis of preheating process 
is valid for all F-term susy hybrid model. For different form of loop correction term just different 
values of parameters have to be chosen and from that exact preheating process can be immediately 
identified from our analysis. 



2 General Formalism of SUSY Hybrid Inflation 

SUSY hybrid inflation models are of two types, F-term inflation and D-term inflation. F-term 
inflation attracts more interest than D-term inflation because it is tailor made to fit with Higgs' 
mechanism. In general F-term inflation assumes a superpotential of the following form[3]. 

W = k$(^-M 2 ) (1) 

Here ^ are chiral scalar multiplets conjugate to each other. $ is a gauge singlet scalar and 
contains inflaton (0). Corresponding F-term scalar potential is as follows: 

V f = k 2 (\^\ 2 -M 2 ) 2 + 2k 2 \<P\ 2 \^\ 2 (2) 

where ip and i/> can be any scalar multiplet conjugate to each other. In the context of applying 
grand unified theories(GUT) ip is thought to be the higgs of the corresponding GUT group. Inflation 
occurs in z/j = ip = trajectory. This is a flat direction with a constant potential of k 2 M 4 . One loop 
correction in this direction gives rise to Coleman- Weinberg potentialpT] of the following form[12j: 

1 M 2 
^) = ^E(-lWlog^- (3) 

Here Fi is the fermion number, M, is the mass eigen value of the component of \1/ and A is the 
renormalization scale. Total effective potential by which inflation is driven, turns out to be 




I (<f> + M ) log — ^ (4) 
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where M is the dimensionality of the representation of ip, ip. In total the potential has two minina, 
one is the local minima at <ft = 0, ip = and another is the global minima at <fi = 0, ip = M. Inflaton 
<f> slowly rolls down through this effective potential during the period of inflation towards the local 
minima. This effective potenial is valid only in the range where <fi > M. After the end of inflation 
another susy-breaking term is required to drive inflaton to <ft = 0. In general an extra mass term in 
the potential, say m 2 2 is assumed for this purpose. But the value of m 2 has to be smaller than the 
co-efficient of <fi 2 arising from the loop-correction k a M 2 . Otherwise this term will effect the values 
of the observed CMB quantities at the time of horizon-exit of a particular mode of inflaton. 
Therefore, in general these types of potentials can be divided into two parts 

V(<P) = V + v(<P). (5) 

where Vo is the constant vacuum energy density and v{4>) is the loop correction term. 

Observed CMB anisotropy parameters are mainly scalar spectral index(n<j) and amplitude of 
density perturbation (A^). These quantities are related to the inflation potential through the 
slowroll parameters(e, rf). 

I - 6e + 2r] (6) 

(T) 
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Here <p in subscript means derivative with respect to <fi and mpi = l/y/G and has a value of 
1.22 x 10 19 GeV. Amplitude of tensor perturbation is related to inflation potential as 



(10) 



k=aH 
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So tensor-to-scalar ratio takes the form r = -A- = 16e. It can be shown that the value of v(4>) is 
negligible compared to Vq, so 
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SUSY hybrid model is a model with negative 77, which gives negligibly small tensor-to-scalar ratio. 
Since no gravity wave signal has yet been observed in WMAP, it puts an upper bound on the 
tensor-to-scalar ratio (r) as 0.36 [13]. If Planck is able to measure any r with a value of 0.1 or above 
this type of model will be discarded. But if it fails to detect r, it will also give an experimental 
restriction on the values of the parameters k, M. 

In standard hybrid inflation scenario <p reaches a critical value <fi = <p c = M, after which waterfall 
occurs and gives rise to the end of inflation. But for this kind of 1-loop potential slowroll parameters 
becomes close to unity before reaches M and inflation ends little earlier than the usual case. We 
have numerically found that a particular relation between k and M will allow to be in the range 
of observed values. 
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3 Different Processes of Preheating 



Mechanism of preheating is a momentum dependent process, i.e. for inhomogeneous field different 
Fourier component evolves differently with time. In parametric resonance processJH] for some values 
of the wavenumber resonance occurs, while for others values it may not. In tachyonic preheating 
[T5] different modes experience different force term in the equation of motion depending on their 
wavenumber. In general both of them gives solution depending on wavenumber of the Fourier mode. 
But tachyonic preheating if encountered by a Fourier mode is lot more effective than parametric 
one. 



Parametric Resonance 

In our particular model inflation ends much earlier than <p reaches the critical value M. So for 
certain ranges of values of M and k inflaton will sustain an oscillation around local minima <p = 
0,^ = 0. Now we Fourier decompose the spatial part of ip as Y^k ipk{t)e lkx where k is the co-moving 
wavenumber. As usual inflaton <fi is regarded as spatially homogeneous. We derive the equations of 
motion for these two fields from equation(T5]) and taking mean-field approximation we get 

k 2 

$ + 3Hij; k + — ij; k + 4« 2 <^ 2 >Vfe + 4« 2 0Vfe - 4k 2 M 2 ^ = (12) 
cr 

4) + 3H0 + 2n 2 (4> 2 )<f) + m 2 <j) = (13) 
(ip 2 ) is the spatially averaged value of ip 2 . By taking X k = a 3 ^ 2 ip k we re-write equation ( II 2p as: 

AX 2 ) 



x k + 



r.1 



2i2 



Ak' 



Ak 2 M 2 



x k = o 



(14) 



Here we have neglected the pressure term since at the time of free oscillation it is vanishingly small. 
Frequency of oscillation of the <p field, oj^, can be estimated from ([TBI at a particular moment as 



^4> — 
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(15) 



where cu 2 = 2k 2 (i(j 2 ) + m 2 . For some ranges of M, k values of is real and for some ranges it will 
be imaginary. This means that for the real values of we can get a sustainable oscillation in </> 
field and for imaginary values the oscillation will be over-damped, i.e. <fi will quickly settle at the 
local minima. 

For oscillatory case of <fi we assume for a certain period = 0o cos u^t and change the variable 
t to r = Uffft. This way we arrive at the familiar form of Mathieu equation[5] from (fT2~|) 



d 2 X k 
dr 2 



+ [A k + 2q cos(2r)]X k = 



where 



k 2 /a 2 + 2k 2 <g + - Ak 2 M 2 



(16) 
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Solution of Mathieu equation takes the form like Xk ~ e UkT where v\, is the critical exponent 
and function of Ak and q. In an instability band of Mathieu chart v gets a real value which gives 
exponential growth in Xk- Occupation number of ipk can be calculated from the solution of eqn. ()16p 
or eqn. ([T2"j) . Since eqn. ffTB]) assumes a purely oscillatory solution of (p we will use numerical solution 
of eqn. f[T2l to get the occupation number from the well known formula 

(19) 

where lo\ = ^ — 4/t 2 M 2 and nk is proportional to \X%\ ~ e 2ukT . 

In an expanding universe parametric resonance is always expected to be in broad resonance 
regime. Values of Ak and q both change with time and parametric resonance is not constrained in a 
single instability band of Mathieu chart. So by observing the growth of n k we can find the average 
critical exponent. In the A, q parameter space broad parametric resonance region is identified to be 
in the area under A < 2q + [7] . This part of the parameter space can give efficient parametric 
resonance. 




Tachyonic Preheating 

Tachyonic preheating is a process in which if a field has tachyonic instability around some extrema, 
its amplitude can increase explosively and preheating gets complete with a single oscillation of 
that field or even before a single oscillation. In our particular model as the inflaton <p crosses the 
critical value M, potential for ip becomes tachyonic. Amplitude of ipk keeps increasing as long as 
k 2 +4K 2 ipl + 4:K 2 (f) 2 — 4k 2 M 2 remains negative. Spontaneous symmetry breaking gets completed after 
it reaches M. Occupation number of the ipk can be defined as equation ffT9~]) . Since col = ^ ~ 4« 2 M 2 
is negative for ^ < Ak 2 M 2 we use uj\ = ^ without any loss of generality [8]. If initially ipk and 
(ip 2 ) is small compared <p and M, we can have 

ipk ~ exp(t v /4«: 2 M 2 - k 2 - Ak 2 ^ 2 ) (20) 

Damping term is neglected for the time being. <j) decreases with time and (tp 2 ) increases very fast. 
So, within a short time growth of ipk and (ip 2 ) gets complete. In the same way initial growth of 
occupation number can also be estimated to be 

n k ~ \tp 2 k \ ~ exp(2t v /4K 2 M 2 - k 2 - 4k 2 <P 2 ) (21) 

So far tachyonic preheating was analyzed neglecting the oscillating cp field. Here we take a 
different approach to analyze tachyonic preheating and parametric resonance in a single frame. We 
use initial form of (ip 2 ) as 

(1P 2 ) = J J ^exp(2t v /4/€ 2 M 2 - k 2 - 4k 2 <P 2 ) (22) 

= — *— ( i + {2uUk 2 <p 2 - Ak 2 m 2 - iyW^-^MA (23) 

Then after some time when (ip 2 ) reaches M 2 we assume that it settles down at that value. 

Since Hubble parameter (H(t)) is also a function of the fields and their derivatives, this as- 
sumption allows us to see the effect of expansion of the universe at the time of preheating. Many 
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authors used to assume static universe at the time of preheating to get read of this complicacy. 
Hubble parameter acts as a damping term in the equation of motions whose effect increases as the 
parameter M and k increases in the potential. Physical wavenumber changes with time as 



k 



k 



(24) 
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At the time of the start of preheating a has been taken to be 1. We did not assume any particular 
distribution of the initial amplitude of the modes. For all modes, initial amplitude is taken to be 
very small compared to <fi and M. Our analysis can be done with any form of initial amplitude 
distribution. Although we find that initial amplitude cannot change the process of preheating; it 
can only change the time taken by the particular mode to start growing. 

The left panel of Fig. [T] and Fig. [2] shows the plot of different quantities for M = 1.22 x 10 15 GeV 
and k = 0.001 and right panel corresponds to M = 1.22 x 10 GeV and k = 0.01. The first row 
in Fig. [1] shows how <f>(t) and (ip 2 ) behave with time. Next we solve the coupled equat ions ( 1T2"TT3"|) 
for Fourier component ipk for k 2 = 0.1(4k 2 M 2 ) and k 2 = 0.3(4k 2 M 2 ) using the above form of (ip 2 )- 
We find both the modes to oscillate around ipk = (not shown in the figure). This goes contrary to 
the standard expectation of parametric resonance that ipk will oscillate around the global minima 
ipk — M. The reason behind it is the oscillating part in the potential coming from the oscillation 
of 4>. We see that among these two parameter sets for the first one damping effect is not so visible 
but for the second set of parameters </> decays rapidly due to the effect of high H(t). 

In second row of Fig. Q] we plot H(t) with time. This H(t) is used in equations (I12|ll3p . We 
see due to the damping effect of the field <f>(t) effects the H(t) decreases with time rapidly for the 
second set of parameter, whereas for the first set of parameter average value of H (t) remains almost 
constant. 

In Fig. [2] we plotted the growth of occupation number(rifc) with time for two different co-moving 
wavenumbers. From the features of log(rifc) we can easily identify two different processes. The lower 
k mode in the first set of parameters gets a steep increase initially. After the start of the oscillation 
of 4> tachyonic preheating ends but parametric resonance becomes effective. For the lower k mode in 
the same parameter set we do not see any significant tachyonic growth in nk, rather the parametric 
resonance becomes effective from start of preheating. In the second parameter set we see that 
tachyonic preheating is effective in both k modes. There is almost no effect of parametric resonance 
due to the high damping effect of large H(t). k 2 = 0.3(4k 2 M 2 ) should not have undergone tachyonic 
preheating because initially it does not experience negative curvature term in the potential. But the 
due to large H(t) physical wavenumber decreases very fast which gives rise to negative curvature 
term as well as tachyonic growth for this mode also. 

We can analytically distinguish between the two processes from the log(nfc) plot of both panels 
in Fig. [2J In the left panel we look at the Mt range between to 20000 for k = 0.1(4k 2 M 2 ). We 
know the amplitude of (ip 2 ) is negligible compared to that of <fi in this range. So curvature of the 
potential remains negative in this range for this particular mode. Therefore ipk is in the form of 
equation fTSOj) . So the slope of log(nfc) with respect to Mt should come as 1.68k. From the plot we 
get this slope to be 1.25 x 10~ 3 in this Mt range. So this result is close to the expected value. 
This assures that the process is tachyonic in nature. This tachyonic amplification should not be 
confused with the narrow parametric resonance shown for some model in Ref. [5]. Now we look at 
k = 0.3(4k 2 M 2 ) line of the same plot for Mt range to 30000. We find that curvature term for 
this mode is positive in this range. This gives an oscillating solution of ipk- So the gradual growth 
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Figure 1: <j)(t) (solid line) and (ip 2 ) (dashed line) are plotted in first row. Hubble parameter H(t) 
is plotted in second row. Values of the parameters are k = 0.001 and M = 1.22 x 10 15 GeV (left 
panel) and k = 0.01 and M = 1.22 x 10 17 GeV (right panel) 
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Figure 2: Growth of occupation number nk for k 2 = 0.1(4k 2 M 2 ) (solid line) and k 2 = 0.3(4k 2 M 2 ) 
(dashed line) has been plotted with k = 0.001 and M = 1.22 x 10 15 GeV (left panel) and k = 0.01 
and M = 1.22 x 10 17 GeV (right panel) 

of log(rifc) is due to parametric resonance for this k. Therefore for these values of M and k both 
the processes are present in different modes. 

In right panel of Fig. [2] we do the same analysis and we find that for both values of k tachyonic 
preheating is the main mechanism. Due to high damping factor the amplitude of ip). gets so much 
damped that parametric resonance can not become effective. We know that the difference between 
log(rifc) for different values of k can give rise to many interesting results like the parameters of non- 
Gaussianity and nonlinearitv|16[ [T7] H~8] . In this analysis we see that this difference gets changed 
depending upon the parameters of the theory. Where in the left panel of Fig. |2] average slope of 
log(rifc) is 1.66 x 10~ 3 and 0.5 x 10~ 3 for low and high k at intial time, in right panel it is 0.007 and 
0.004 for the corresponding k values. 

There is a distinct difference between the time scale of the right panel and left panel. This 
comes from the frequency of the oscillation of field in equation ([15]). Since Uq is proportional 
to k 2 and H 2 is also proportional to k 2 time period varies inversely with k. This observation is also 
in accordance with ref.[8], where time taken by the field to reach the zero curvature point in the 
potential was expected to vary as 1/ n 

For analyzing preheating in SUSY F-term hybrid model potential of the following form was 
assumed [6] 

V(i/j) = k 2 M 2 + 3/«V - 4k 2 i/j 3 M (25) 

Potential in eqn.(|5]) takes this form if we assume a trajectory for fields as ip = ±(0 — M). Growth 
of the perturbations with time had been shown for this potential using LATTICEEASYjH]. But 
different mechanisms are not unidentifiable from this work. In our analysis two different mechanisms 
are clearly identifiable. 

4 Results and Discussion 

We have shown in Fig. [3] that there is relation between k and M for the particular observed value 
of which is (2.41 ± 0.11) x 10~ 9 . To do the analysis we have taken the number of e-foldings, N, 
corresponding to the WMAP pivot scale 0.002 MPc" 1 , to be equal to 60. And M has been taken to 
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Figure 3: The curved line in k, M space belongs to the values of k and M for which A 2 ^ = 
(2.41 ± 0.11) x 10~ 9 . The blue part of the line corresponds to the n s = 0.983, the lowest value 
possible for the model. 

be 8. For a particular value of M there are two possible values of k. So to choose the appropriate 
value of k and M we have to look at the other observables. Points in this curve correspond to 
different values of n s and r. For all the points, the value of r is of the order of 10~ 7 , which is 
negligibly small as expected in negative rj models. n s achieves its lowest value of 0.983 around 
M = 9.55 x 10 15 GeV and « = 0.008 - 0.01. This is a problem with the standard SUSY-hybrid 
model. This lowest value of n s is not in accordance with the latest measurement by WMAP, which 
gives n s = 0.963±o*^2- To fit with this observation various correction terms are added in this effective 
potential [IHl [20]. In all those models for the observed value of A^ a relation between k and M can 
be found. As we have said earlier we chose the standard form of the susy-hybrid inflation just for 
an example. 

We have taken this tuned values of parameters in our analysis of preheating and the results are 
plotted in Fig. gj We see that both the k 2 = 0.1(4k 2 M 2 ) and k 2 = 0.3(4k 2 M 2 ) behaves almost in 
the same fashion of k = 0.001 and M = 1.22 x 10 15 GeV case. Lower k mode undergoes tachyonic 
preheating initially and higher one undergoes parametric resonance. Average of the critical exponent 
( z/ fc)avg for the higher k mode has been estimated from log(nfe) plot of this figure and it comes to 
be 0.22. For the lower k mode growth rate of log(rife) is 0.012 with respect to Mt. This parameter 
fixing is necessary for the cases where gravitational wave production is predicted from preheating 

[25- 

Broad parametric resonance is generally thought to occur for the oscillation of <p around global 
minima. In that case value of ip is around M and amplitude of <fi is also of the order of M. 
Here we have seen that just after the end of inflation if inflaton can sustain oscillation around the 
global minima preheating via broad parametric resonance is possible. For higher values of k and 
M amplitude of (f> and iftk decays down quickly which helps tachyonic preheating to dominate over 
parametric resonance. 

Standard susy-hybrid model has only two parameters k and M. But other models may involve 
some more parameter. In that case it might be hard to pinpoint the values of parameters for which 
it will give two different preheating process. But this kind of models has been highly popular in 
the context of using Grand Unified Theories (GUT) in inflation. It was expected to get M as the 
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Figure 4: In left panel (j)(t) (solid line) and (ip 2 ) (dashed line) are plotted for k = 0.01 and M = 
9.55 x 10 15 GeV. In right panel, for the same values of parameter, growth of occupation number 
has been plotted for k 2 = 0.1(4/t 2 M 2 ) (solid line) and k 2 = 0.3(4k 2 M 2 ) (dashed line). 

vacuum expectation value of higgs multiplet from the symmetry breaking scale of a particular GUT 
group, ip will have the dimensionality of higgs. If those kinds of attempt are successful in future 
it will restrict the value of M. So, in those cases even with higher number of parameters it will be 
possible to tell about the exact preheating process. 

Exact dynamics of preheating is necessary to understand for studying nonlinear evolution of 
density perturbation [22] . Imprint of preheating can be observed in density perturbation of large- 
scale structure which originates from the super-horizon modes of fluctuation of the scalar fields at 
the end of inflation[23j. So, to find out the best possible model of inflation, along with reproducing 
other CMB observables, analysis of preheating process is also important, first step to which can be 
identifying the actual process. 

Acknowledgments: We thank Palash B Pal and Koushik Dutta for discussions and valuable 
suggestions. 
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